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Abstract 

The variety of nanowire crystal structures gave rise to unique and novel transport phenomena. 
In particular, we have explored the superlattice profile generated by strain field modulation in twin- 
plane nanowires for the tuning of transport channels and the built-in spin-orbit potential profile of 
polytypic nanowires, in order to realize a spin filter. The Multicomponent Scattering Approach has 
been used in terms of the Transfer Matrix Method to describe the phase-time of charge carriers. 
This system showed advantages for attaining conditions for the propagation of wave packets with 
negative group velocity. Moreover, the spin transport effect of a potential profile with volumetric 
spin-orbit bulk inversion asymmetry, as present on polytypic nanowires, was described through 
the Reverse Runge-Kutta Method. Using the peculiar symmetry of the excited states we have 
characterized a dominant spin dependence on structural parameters that results in effective spin 
filtering. 
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I. INTRODUCTION 


The integration of nanowire (NW) architecture as a building-block for electronic and 
photonic applications has become an issue of considerable research interest.-*^ Therefore, the 
study of their transport properties is a central task for both theoretical and experimental 
endeavors boosted by the developments of growth techniques that allow minute electronic 
structure engineering. The ability of thorough control of NW crystal structures by varying 
parameters as growth temperature and pressure, NW diameter and surface density, precursor 
molar fraction, III/V aspect ratio, and incorporation of impurity atoms allows the fabrication 
of a variety of high quality wires displaying zinc-blende, twin-planes (Fig. [11(a)), stacking 
faults (Fig.[T](d)), and wurtzite crystalline structures.-^- The improving of this level of control 
can lead to a variety of electronic behaviors and these new prospects have opened different 
opportunities for the tuning of their transport properties. 

The effect of volumetric spin-orbit (SO) bulk inversion asymmetry (BIA) in wires with 
zincblende crystalline structure has been previously discussed,- where the degeneracy of 
spin states for each Bloch wavevector k, in the absence of an inversion center, is broken. 
This leads to the spin splitting of energy states of both electrons and holes without lifting 
the Kramer degeneracy in the absence of an external magnetic held. It was shown that, 
although the ground state of a cylindrical NW is always doubly degenerate in terms of the 
spin polarization at any value of k^, regardless of the direction, different terms induced by 
the Dresselhaus SO contribution, split excited levels when considering terms proportional 
to kz, the longitudinal component of the wavevector. This is a dominant process that would 
preserve the spin polarization for small values of kz opening up the possibility of using 
transport measurements to explore the preferred spin channels according to the direction of 
propagation of spin polarized currents along the quasi-ID nanostructures.- In this way and 
by using the hexibility of manufacturing polytypic nanowire (PNW) structures, we developed 
transport simulations in a NW presenting regions with and without centers of bulk inversion 
asymmetry that act as spin scatterers. Thus, if carries are injected along the z-axis of a 
PNW with certain incident energy, the peculiar symmetry of excited states in the zincblende 
cylindrical NW region allows only the transmission of a high degree of spin-polarized current 
according to the propagation direction. A dominant spin-up or spin-down character of the 
current depends on the sign of kz and the current density is directly linked to the interplay 
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FIG. 1. (a) Schematic diagram of a polytypic nanowire (PNW). (b) Electronic structure profile 
of a PNW used in the simulation of the spin transport properties, where /3(z) is the Dresselhaus 
or BIA spin-orbit potentials, (c) Schematic diagram of a twin-plane nanowire (TPNW) used in 
our simulation for the electronic structure, (d) Electronic structure prohle of the top-most valence 
subbands, (e) Notation used to characterize the scattering potential profile. The segment [zL,zji] 
represents the periodic cell^ for the traveling modes. 

between the wire radius and the wavevector. 

Likewise, motivated by the successful conhrmation of surface effects on transparency 
modulation of ID twin-plane based superlattices (TWPSL), the Multicomponent Scatter¬ 
ing Approach (MSA)^>^ was applied to model and prescribe appealing tunneling events for 
both uncoupled light-(Ih) and heavy-holes (hh) through these NWs. Phase-time calculations 
show several anomalous transport properties for each kind of charge carriers and reflect, as 
expected, the mini-band spectrum for ID-TWPSL.— For instance, a clear hallmark of neg¬ 
ative values of the phase time will be described. Undoubtedly, such evidences are known 
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and were reported before for electrons, holes, and twin-photons in typical III — V layered 
structuresii^— but, as far as we know, they are novel events for TWPSL wires. Addition¬ 
ally, under certain circumstances, oscillating regions were detected for both carriers as they 
trespass the n barriers of the ID-TWPSLs. This slightly recalls the Ramsauer-Townsend 
effect 

By heeding these two guidelines, the rest of the paper is organized as follows. The 
next two sections present the theoretical methods used to tackle the transport simulations 
through spin-orbit interaction centers and twin-planes lattices. Afterwards, the main results 
are presented followed by the concluding remarks. 


II. DRESSELHAUS SPIN-ORBIT EFFECT IN POLYTYPIC NANOWIRES 


We simulate various transport properties in a NW presenting regions with and with¬ 
out centers of inversion asymmetry (BIA) introduced by the Dresselhaus SO potential as 
sketched in Fig. [I](a). The wavevector has to be considered as an operator which fulfills the 
temporal inversion symmetry and the single-particle Hamiltonian for the conduction band 
of a NW, including Dresselhaus SO interaction (SOI), can be written as H = Hq + "Hd, 
where Ho = /2m* + V{r,z)]I is the Hamiltonian for uncoupled spin-up and spin-down 

states, I is the 2x2 identity matrix, and V (r, z) is the total spatial confinement with lateral 
r = {x, y) and longitudinal z along the wire axis for the set of cartesian coordinates. The 
effective BIA Hamiltonian can be separated into three contributions. 


Hv = 



( 1 ) 


where Hiv = Ii 2 v = > and H^v = |^+, (^+- 

They are, respectively, the linear, the quadratic, and the cubic Dresselhaus SOI contributions 
written in terms of operators k± multiplied by the dimensionless parameter 7 = lD/{o?£-a)-, 
and with 80 = h^/{2m*a^) being the energy scale for the NW conhnement inside the cylin¬ 
drical region. The Schrodinger equation for H is separable in the entire level spectrum as a 
product of purely in-plane localized function of the radial coordinate r times the free motion 
function along the 2 ;—axis, and can be written as. 
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/ (A-+ih,^P{z)-§; \ / fii(.z)\ ^ / -B 0 ^ / f"A^) 

V -j-^Az) + 9tMz)& (A. - gr)^„. - ih,yD(z)§; ) f‘J^(z) ) V 0 B j \ /^(z) 


Here, the index s represents the set of quantum numbers (n, m), ordered by increasing values 
of energy with ^ [Eq/Eq — a?k1) being the m*^-zero of the n^^-order Bessel 

function, {fin,m) = 0, and with the index a accounting for the spin-up and spin-down 


states. The innermost matrix elements are gf , = (s 


± 


^ / 7 ^ss' 2 ^ ^ 


kl + kl 


s'), and 




{klAkl-kl} 


The Dresselhaus coefficient, as function of position along 


the wire axis, can be emulated as 


+ -erf(^-^)], (3) 

where erf(z) is the error function and 2 £ is the width of the resulting quasi square region 
with hnite hrst derivative, illustrated in Fig. [T] (b). Note that the method presented here 
can afford any kind of potential regardless its spatial shape. After symmetrization of the 
Dresselhaus Hamiltonian, in the dimensionless form, one reaches a nonlinear set of equations, 
written as 


- fa/A) + iPiAhfaa'A) + ^(^9 

- 5 {inz)hfAjz) + {P"(z)g- + 2j+/3(j)) /•f;(z)) 


and 


- + mrc-'iz) 

+A{z)rfA^{z) = {e+^A)rJAz) 

+ i {iP'{z)hf^i{z) - (/ 3 "(z)r - ^J-Piz)) f^::{z)) , 


( 4 ) 


( 5 ) 


where k = ak^, = agf^,, h = a^hss' and order to solve this set of equations 

we have used the Reverse Runge-Kutta method as shortly described below. 
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A. Reverse Runge-Kutta 


For two spin channels and a single subband s = s', one can write the wavefunction in 
each of the three regions delimited in Fig. [T] (b) as, 

^—izk JD I ^izk 

e Kaa + 6 

-.—izk ; 


d'T = 


e 


^II = 


, and \l/iii = 


^izk'Y' 

^ ^ OLOL 

^izk'~r 


respectively. From the right interface boundary condition, \l/ii(F) = \l/in(F), one can write 

Azk 


Fiii(L) = 


Fii(^) = 


1 AoiOii^L') \ 

/ 

I Taa \ 

- 

\ 1 

■ 

V '^a’a ) 

V 

1 faaiy ^ 


\ faUy , 



Azk 


, and 


1paa{z) 

Tea 


where Fj is the dispersion region function. Hence, Fn(F) = Fni(F), and Fjj(L) = Fjjj(L) 
are, respectively, 

ALk 


1 . I L"(L) 

, and 

,<Lt I 1 


iLkj^ 


%e 


( 6 ) 


mL) 
faAL) 

Following the reversal procedure for the left boundary condition, where \l^i(—L) = 
\l/n(— L)), one can address the unknown function Fi as, 

'^aoi ( -^) 

= I I ■ 


Fi(2) = 


Tolol 


Therefore, from the left boundary condition, Fi(—L) = Fn(—T), and Fj(—L) = Fjj(—L), 
along with the current conservation, |Tq,q,(A:)|^ + \Raa{.kA + \Ta'a{kA + \Ra'a{.kA = 1, one 
is able to dehne the transmission (T) and reflection (i?) coefficients in terms of the spin 
components. Thus, 


Taa{k) = 


2 k 


kf^s^{-L)-tfjA-Ly 


(7) 
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( 10 ) 


( 8 ) 


(9) 


where the f{—L) is obtained employing the right boundary conditions, Eq. ([6]), along with 
the 2 ;-component Schrodinger Eqs. (jl]) and ([5]). Similar expressions are obtained extending 
this method for multiple subbands. 

III. PHASE-TIME IN TWIN-PLANE SUPERLATTICES 

In turn, for the simulation of the tunneling properties in the twin plane NW we have 
used the transfer matrix framework assuming a ID-TWSL, represented in Fig. [T](c) with 
alternated layers whose quantum heterogeneity, due to strain held modulation, is revealed 
along the 2 ;-axis in Fig.[T](d). The envelope function (EE) coefficients at both TWSL extremal 
slabs, are bounded as 


A 


Mi(Zo)-'-M2(z2)...-M4(z4)”'- 


( 11 ) 




M5(z4) 

B, 


being 



the transfer matrix. The EE is given by 



( 12 ) 


and its weighted derivative, as 



( 13 ) 
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at both layers / and V of the ID-TWPSL. [see Fig.[T][d)]. The boundary conditions for scat¬ 
tering quantities can be set by assuming an impinging mixing-free particle (quasi-particle) 
traveling from the left layer I, represented in Fig. [I][d), 


Ai — 1 ; Bi — T 
A 5 = R ; i?5 = 0 


(14) 


where T and R stand for transmission and reflection amplitudes, respectively. Next, we 
establish a correlation between quantities fITT]) and flTT)) . as 


1 

R 


M{zi,Zr) 


T 

0 


I Mn{zi,Zr) Mi2{zi,Zr) 
yM2l{Zi,Zr) M22{Zl,Zr) 



which straightforwardly lead us to 


( 15 ) 


Tij ^ 7 ®., + ^ \Mii{zi,Zr)] 


-1 


-o 


Rii — R^ii + iR^zz — 


M2l{zi,Zr) 


(16) 


Mxx{zi,Zr) ’ 

the amplitude matrices for transmission and reflection, respectively. Next, we obtain from 


IT- 


^3\ 


and 


the scattering probabilities for the transmitted and reflected fluxes, re¬ 
spectively. For the calculation of the tunneling phase-time, a simple procedure is followed. 
Firstly, we calculated the transmission phase 


6 ij = arctan 


Tuz 


+ f (1 


\Tuij I 


( 17 ) 


arctan 


\»(Mll(2i,2,)) 


+ 


1 - 


5)?(Mll(2i,2r)) 

|K(Mll(2i,2,))| 


and afterwards, we derive the n-cell tunneling phase-time^ii^ (pjOn “ assuming left¬ 

incoming particle (quasi-particle) 


(Ti)n 





(18) 


IV. RESULTS AND DISCUSSION 


The formalism developed in Sec. |TT] allows the study of spin transport of PNWs, where 
changes in the zincblende region width, w, the NW diameter, 0, the incident energy, T, and 
















the SOI strength, /9, can be used as tuning parameters. The presented results correspond 
to an incident spin polarized current injected on the wurtzite region of a InSb PNW, where 
(3 = 760.1 eV-A^. We will see that the SOI within the zincblende region gives rise to a 
small net polarization. Projecting it to several sequential regions, it is possible to realize 
spin hlters and modulate the densities of spin-polarized current flowing in both parallel and 
antiparallel directions of the NW. 

In order to characterize the system response to the injection of a spin-unpolarized super¬ 
position of spin-up and spin-down currents, we calculate the spin persistence ratio as defined 
by [|T„„|2 + )]/T, where T = |T„„|2 + . 

Figure [2] shows the color code maps of characteristic persistence vs. incident energy, £, and 
ta, and vs. £ and (p. Note, that the persistence is always positive, while characterizing the 
transmitted current, but does not reach the unity (Figs. |2] (a) and (c)). This positiveness 
refers to a signature of almost perfect spin-preserving channel, as expected due to the non¬ 
degenerated ground state,- but there are oscillating regions which evolve with changes in 
£, w, and (p. On the other hand, from the characterized reflected current it is possible to 
observe that the persistence goes from spin-preserving to spin-reversing (negative values) 
channels. The spin-reversing regions are identified by negative numbers and also evolve 
with changes in £, w, and p. Despite the difference of spin-preserving and spin-reversing 
channels, the oscillation pattern of the reflection coefficient is doubled in comparison to the 
transmitted one and this should be expected once the current crosses twice the region of 
Dresselhaus SOL— Moreover, the evolution of the oscillations goes with both £ and p, so 
that £p = constant, as shown in panels (c) and (d) of Fig. [2l due to the modulation of the 
energy levels with p, Eq. ([2]). 

When considering excited states in the transport, one can expect a small net spin po¬ 
larization for the current. This is due to the fact that for transport along the ^-direction, 
the Dresselhaus SOI will produces the precession of both spin-components but in different 
ways.- The degree of spin-polarization is defined by [iT^Qp -|- \Ta>a\‘^ — (l^o'a'P + \Taa'\‘^)]/T. 
Figure [3] shows the polarization as function of £ and w in panels (a) and (b), and as a 
function of £ and p, in panels (c) and (d). 

Observe that the transmission displayed in Fig. [3] (a) shows well dehned polarization 
regions which evolve with both £ and w. Around £ = 5 meV and w = 150 A, a smaller 
percentage of spin-down polarization is evident, as well as for = 45 meV. Looking at 
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FIG. 2. Persistence contour maps of polytypic InSb quantum wire for the first band where (3 = 
760.1 eV-A^. Panels (a) and (b) show the transmission and reflection as function of the barrier 
width and energy for a wire with radius 101 A. In panels (c) and (d) we show transmission and 
reflection as function of the wire diameter and energy. 


w = 800 A, one can see an opposite spin-polarization for = 45 meV. Moreover, it is 
clear that at this width ta, a modnlation of the polarization takes place, from spin-down to 
spin-np character, by varying 8 . This last hallmark is also revealed in the reflection panel, 
shwon in Fig.[3](b), bnt with opposite polarization featnres. Additionally, Figs. |3](c) and (d) 
show the same behavior of Figs. |2] (c) and (d), with oscillations occnrring by varying both 8 
and 0, snch that =constant. These results indicate the possibility of modulation of the 
polarization by selecting the incident energy of carriers, the NW wire width, the zincblend 
region width, and the SO potential. 

The presented results for the spin persistence and polarization refer only to a config¬ 
uration of one zincblend region between two wurtzite layers, as seen in Fig. [1] (b), and 
this was enough to characterize a small percentage of spin flux modulation. Increasing the 
number of layers containing the zincblende regions, one can expect an increase of the polar¬ 
ization. Figure m shows the transmission, (a), and reflection, (b), polarizations for multiple 
wurtzite/zincblende/wurzite regions. As can be noted, an increasing number of layers leads 
to increasing spin hltering and enhancement of the hnal polarization. 

Along with the polytypism, the growth conditions may lead to controlled stacking faults 
appearance and the synthesis of twin plane superlattices. For these later, as represented 
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FIG. 3. Polarization contour maps of polytypic InSb NW for the first subband (3 = 760.1 eV-A^. 
Transmission (a) and reflection (b) as function of the barrier width and incident energy for a wire 
with radius of 101 A. Transmission (c) and reflection (d) as function of the wire diameter and 
incident energy. 
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FIG. 4. Polarization of polytypic InSb NW with multiple wurtzite regions, n, where /3 = 
760.1 eV-AA (a) Transmission and (b) reflection as function of the energy for a wire with ra¬ 
dius of 101 A. 

in as Fig. [1(c), we have proposed a transport analysis of the tunneling phase-time for both 
Ih and hh. Here, the results for electrons were not shown due to similarities with the Ih 
transport results.— 

Using the formalism of Sec. IIIIl the mini-band spectrum for holes tunneling through 
(n = 2 - 100) ID-TWPSL is depicted in Fig. 0 as ([IB]) and evolves with the carrier 
incident energy, 8 . As expected, this simulation agrees with an previous direct calculation 
of the ID-TWPSL spectrun>^ and responds to the mini-band spectrum of the ID-TWPSL, 
i.e. it reproduces accurately the quasi-stationary holes states of the embedded ID-TWPSL 
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quantum wells. We found the double resonant tunneling barrier (DBRT) curve {n = 2, red 
dashed line) to represent a lower bound for the ID-TWPSL mini-band spectrum as had been 
reported for electrons^ and for holes^ in two-dimensional systems (see Fig. |S] (a) and (c)). 
Appealing regularity for energy sections has been found, where rioo < t'so < 'r 2 < 0 (see dips 
in Figs.|S](b) and (d)), where the transmission vanishes and, thereby, coincides with the mini¬ 
band forbidden regions of the system (Figs. [5](a) and (c)). In other words, this abnormal 
behavior for holes is directly connected to a large backscattering {\T{S)\‘^ ~ 0).— The 
later could be due to the fact that negative values of r are not only related to the incident 
wave from left but also might be a combination of both incident and accumulatively reflected 
waves from the ID-TWPSL quantum barriers traveling towards the left-hand side.— One 
may notice that the quasi-classical free motion time, Tf = which characterizes the 

temporal scale for a particle with mass rUe// to travel through a none-scattering space region 
of dimension n times the thickness L^,- can be straightforwardly compared with Eq. (fT8|) . 
Indeed, in our case, a striking quality arises whenever {Tf)^ > r 2 , 30 , 100 , for the mini-gap 

regions, respectively (Fig. [5](b) and (d)). Therefore, for incoming hh states, with energy in 
the neighborhood of 20 meV, ('r/)gQQ —tioo ~ 60 ps, while for incoming Ih with incident energy 
~ 120 meV, — Tioo ~ 15 ps. This earlier arrival time for both hh and Ih, previously 

predicted for electrons^S and holes^ that tunnel through a two-dimensional DBRT and a 
SL, suggests a faster passage of holes trough a ID-TWPSL, under specihc conditions. The 
regularity of more speedily passage of hh, compared to that of Ih,- is clearly preserved in our 
numerical simulation. Importantly, only at low energies (~ 5 meV for hh, and [14 - 23] 
meV for Ih) rso approaches (t/) 3 q from upper side down, as detected for electrons^S and holes 
in two-dimensional SL.- Additionally, it seems to be satished that (Tf )2 3 oioo ^ D, 30,100 for 
the rest of impinging energy values in the transparent region, contrary to the two-dimensional 
case.- 

In order to spread some light and try to elucidate some features of the negative phase-time 
(eq. ffT8|) i values, one can see in Fig. [6] the transmission amplitude, = T^.. +zT<^.., in the 
complex plane for different numbers of cells, n. Here, 0ij{S) has been taken as the angle and 
\T{S)\‘^ as the radius. As can be seen, the major signature of the curves is their clockwise 
direction evolution, which can be directly related with positive phase-time values. However, 
counterclockwise segments clearly appear for hh with n = 3 (green-triangles) and n = A 
(blue-diamonds), characterized by loops in Fig[6][b). Interestingly, these loops seems to be 
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FIG. 5. Calculated transmission probability and phase-time {tu) for various sequences of twin- 
plane cells, n, in a ID-TWPSL as a function of incident energy {£): (a) and (b) Transmission 
probability, r and quasi-classical free motion time (r/), respectively, for light-holes (Ih). (c) and 
(d) Transmission probability, r and r/, respectively, for heavy-holes (hh). 

modulated by the n = 1 (black-squares) curve as a lower bound, Fig[H](b). The last could 
be a certain resemblance of a similar average behavior, detected for electrons^^ and holes^ in 
two-dimensional systems. The counterclockwise direction of Oij{S) is directly related with 
negative values of the phase-time. For Ih, in Fig. |6](a), this occurs when \T{S)\‘^ ~ 0 for 
n > 4 and for a large number of cells. One can hnd this negative values of r for hh when 
|T(£l)p ~ 0 for n > 3, in Fig. El^b). Despite some criticism in the specialized literature, 
clear hallmarks of negative values for r are reported, where any numerical artifact has been 
accurately excluded. In spite of the interpretation proposed above, for this counterintuitive 
topic, a conclusive robust theoretical modeling remains a puzzle. Such evidences are novel 
events for TWPSL NWs, although they are well established on other systems.—^— 

Figures [71(b) and (c) display the calculated r^, eq. flTSl) . for different sequences of TWP 
cells caused by the strain effect, as a function of the barrier width, L^,. Oscillating regions are 
detected for both hh and Ih, as they trespass the n = 30,100 ID-TWPSLs with low values 
of energy, £ < 0.514 [see Fig. [7] (b), (c), (e),and (f)]. Accumulative barrier interferences 
with evanescent hole states have been assigned as the possible origin of these oscillations. 
This slightly recalls the Ramsauer-Townsend effect.- This behavior fades away abruptly at 
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FIG. 6. Representation of the transmission amplitude through different numbers of cells, n: Tij — 
Tjf , in the complex plane for increasing incident energy, (a) Ih and (b) hh. Note that both 

diagrams have been rotated for clarity. 


Lft 50 A for Ih whereas for hh is observed at Lj, 30 A. It can be argued that beyond a 
critical value of barrier width, the wavefunction penetration length becomes negligibly small 
and, thereby, the \T{£)\‘^ tends to vanish while < 0 values arise. None of such oscillations 
has been found for the DBRT case (n = 2). 



FIG. 7. Calculated for various sequences of twin-plane cells in a ID-TWPSL as a function of 
the barrier width (Fb). For Ih: (a), (b) , and (c). For hh: (d), (e), and (f). 
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V. CONCLUSION 


In summary, motivated by well established methods for manufacturing polytypic and 
twin-plane NWs, we have proposed two different carrier transport studies for the charac¬ 
terization of these systems. The hrst one was the spin transport simulation on a polytypic 
NW structure presenting regions with (zincblende) and without (wurtzite) Dresselhaus SO 
potential. For this transport calculation we developed a new approach based on the Reverse 
Runge-Kutta method. We have found that carries injected with certain energy along the z- 
axis of a TPNW, would allow the transmission of spin-up current along the -fz-direction and 
spin-down current along the — 2 ;-direction, as due to the peculiar symmetry of the excited 
states in the zincblende cylindrical region. Furthermore, the increasing number of layers 
in the structure leads to increasing spin hltering ability as well as an enhancement of the 
hnal current polarization. The spin-polarized current density can be directly linked to the 
interplay between the wire radius and width of the Dresselhaus SOI layer. 

In turn, the second simulation uses the MSA and TMM methods to study the charge car¬ 
riers passage through a strained ID-TWPSL. As expected, the calculation averaged out the 
mini-band spectrum for ID-TWPSL. Amplitude oscillating regions were obtained for both 
hh and Ih carrier currents with low values of incident energy, thus evoking the Ramsauer- 
Townsend effect. Moreover, the phase-time assumes negative values whenever the ID- 
TWPSL properties reach a minimum or negligible transmission probabilities. This abnormal 
behavior is directly connected to a large carrier backscattering at barrier regions of the hnite 
superlattice. 

In order to elucidate this effect, we have characterized the transmission amplitude in the 
complex plane where the major clockwise signatures were directly related to positive values 
of Tn, whereas the minor counterclockwise ones were connected to negative values of the 
phase-time. This effect is only observed for structures with n > 2. We hope that these 
hndings would stimulate further search for the realization of spin hltering structures. 
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